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ABSTRACT 

The radiative efficiency of AGN is commonly estimated based on the total mass accreted and the 
total AGN light emitted per unit volume in the universe integrated over time (the Soltan argument). 
In individual AGN, thin accretion disk model spectral fits can be used to deduce the absolute accretion 
rate M, if the black hole mass M is known. The radiative efficiency 77 is then set by the ratio of the 
bolometric luminosity Lboi to Mc^. We apply this method to determine 77 in a sample of 80 PG 
quasars with well determined Lboi, where M is set by thin accretion disk model fits to the optical 
luminosity density, and the M determination based on the bulge stellar velocity dispersion (13 objects) 
or the broad line region (BLR). For the BLR-based masses, we derive a mean log 77 — —1.05 ± 0.52 
consistent with the Soltan argument based estimates. We find a strong correlation of 77 with M, rising 
from 77 ~ 0.03 at M = IQHIq and L/L^^d ~ 1 to 77 - 0.4 at M = IQ'^Mq and L/iEdd ~ 0.3. 
This trend is related to the overall uniformity of iopt/iboi in our sample, particularly the lack of 
the expected increase in iopt/^boi with increasing M (and decreasing L/LEdd), which is a generic 
property of thermal disk emission at fixed r/. The significant uncertainty in the M determination is 
not large enough to remove the correlation. The rising 77 with M may imply a rise in the black hole 
spin with M, as proposed based on other indirect arguments. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — galaxies: 
quasars: general 



1. INTRODUCTION 

Material falling in nearly circular orbits onto a black 
hole (hereafter, BH) looses a fraction of its rest mass en- 
ergy during the infall. The lost energy is emitted as an 
outflow of radiation and particles (and potentially Poynt- 
ing flux) . A measurement of the fraction of mass inflow 
M converted to radiation iboi, provides a measure of the 
radiation efficiency 77 = L-bo\/Mc^ . 

In the "standard" ac cr etion disk (AD) model 
(IShakura fc SunvaevI ITQTI INovikov fc Thornd IT973I) . 
the BH spin a* determines 77 because it sets the 
marginally stable orbit, rms, beyond which the ma- 
terial is assumed to fall into the BH without loosing 
further energy. Since the total efficiency is a rising 
monotonic function of a,, the measured 77 provides 
a lower limit on a*. The value of a* is important 
because it tells us how the BH mass M grew. If it 
grew mostly through a single event (major merger 
or continuous gas accretion) then a, will be close to 
unity. If it grew through a series of independent events 
(minor mergers, episodic accretions), then a* will be 
close t o zero iH ughcs & Blandford 2003; Gammi ^ et al.l 
2OOI IVolonteri et al.. .2005; .King fc Pringla 120061: 



Berti fc Volonterill200g 



In the absence of torques near rms, the value of 
a* sets, through the value of Tms, the spectral en- 
ergy di stribution of the accre ti on di sk (iCunningham 
T975I IKolvkhalov fc SunvaevI HOSl ISun fc Malkan 



Tosa ILaorfcNetzeri 1198a ISinceh fc KrolikI 119981: 
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IHubenv et al.|[2000l) . the rotation of the polarization an- 
gle of the A D emission, (Connors et al.' 1980: Laor et alj 
Il990t iDwBiak et al.. ,2008: Schnittman fc Krolik 2009i) 
and the profile of lines emitted by the AD (iFabian et al.l 
19891: IK oiimallTggTI: iL aoIl [T99lt 'iDabrowski et al.' 119971: 
Beckwit h & Done 20041: IBrenneman fc Revnol ds 2001 
Revnolds fc Fab ian 2008^). These methods are currently 



limited by the available quality of the data, and by po- 
tential uncertainties in our models of the AD structure. 
As a result, we do not yet have precision measurements 
of a, in more than a few objects. Thus, an additional 
independent constraint on a* based on ?7, is useful. 

A determination of r/ can potentially provide an up- 
per limit on the additional power which may be gen- 
erated by the accretion in a jet/wind outflow. With- 
out torques, there is an upper limit on the total ef- 
flciency of 40% for a, — 1, or 31%, for the maxi- 
mal spin within an AD of a, = 0.998 (jThornd [1974). 
Such outflows are important as they can couple to 
the surrounding gas more efficiently than radiation, 
and may significantly affect the host galaxy evolution 
(e.g. IMcNamara fc NulsenI 120071), suppress cluster coo l- 
ing flows re.g. IChurazov et arn200l lAllen et al.l [2001 . 
and may be relevant to the correlati on of the black 
hole mass with the galaxy properties (jMagorrian et al.l 
119981 and citations thereafter). The implied jet power 
of AGN in cooling flow clusters can be significantly 
larger or smaller than the radia tive power, depend- 
ing on the AG N luminosit v fe.g. iSh ankar et al.l 120081 : 
iMerloni fc Hehi z 2008; Catt aneo fc Be st 2009). Clearly, 
it is useful to get an independent upper limit on the 
ratio of mechanical/radiative power, based on a direct 
determination of 77. If magnetohydrodynamics torques 



(IGammielll99aiKroliklll999l:iMcKinnev fc Gammidl200l 
iDe Villiers et al.l l2003| ) are present, then the maxi- 
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mum efficiency can (instantaneously) exceed the limits 
for a no-torque dis k and even exceed unity (see e.g. 
lAgol fc Krolik|[2000D as the flow taps the spin energy of 
the BH. Therefore, credible estimates of such large ef- 
ficiencies would provide evidence that such torques are 
p resent i n real accretion flows. 

iSoltanI H1982D noted that the global AGN average ra- 
diative efficiency, ryav, can be estimated for the AGN 
population by comparing the integrated M per unit vol- 
ume at the current epoch, with the in tegrated AGN 
luminosity per unit volume over time. I.Soltan! (Il982f) 
also showed that ryav is elegantly independent of the 
cosmological model (a major unknown at that time). 
Recent studies based on the Soltan argument lead to 
>v > 0.1 (e.g. lYu fc T remaine"2002': 'E lvis et all [20021: 
IMarconi et all 120041: IBargcr ct al. 2005). This method 
has also been used to estimate the time and lumi- 
nosity depen dence of 7?av t hrough more detailed mod- 
eling (e.g. iHopkins et all 2006; Shan kar et all 120091: 
I Wang et al.ll2009l:lRaimundo k Fabian..2009l ). but the de- 
rived values are significantly uncertain. 

The purpose of this paper is to discuss a method to 
derive rj directly in individual AGN. The method as- 
sumes that the optically emitting regions of QSOs are 
accretion powered and radiatively efficient, thus grav- 
itational binding energy is dissipated and radiated lo- 
cally within the AD. The corresponding thin AD mod- 
els were calculated to increasing levels of details, from 
the simple local blackbody approximation to stellar at- 
mosphere like models where the vertical structure and 
the local s pectrum are calcula ted with increasing ac- 
curacy fsee iHubenv et al] 120001 and references therein). 
The integrated thin AD luminosity density turns out 
to be largely set by M and Al . Thus, one can derive 
M based on the observed Ly, if M is known. This 
met hod has been u sed p reviously by ICoUin et al.l H2002f ) 
and IBian fc Zhaol (f2003l) . using simple analytic expres- 
sions, valid at long wavelengths for the emission of an 
Newtonian, thin, blackbody AD (e.g. iBechtold et al.l 
to determine M in a sample of AGN. iCollin et al.l 
assumed a value of rj to estimate Lboi and in- 
ferred that many low M AGN must be super-Eddington 
accretors. Bian & Zhac^ ()2003D estimated Lboi indepen- 
dently for each object, which they then used to estimate 
77, yielding an average log 77 — —1.77 ±0.49 in a sample of 
radio-quiet AGN, and log r] = —0.90 ± 0.62 in radio-loud 
AGN. 

Observations indicate that simple thin AD model can- 
not reproduce the overall SED. This is due to reprocess- 
ing (IR), Comptonization in a corona (X-ray), radiative 
transfer effects in the inner AD, a thick AD, etc'. Our 
method relies on the viability of the simple thin disk ap- 
proximation in the relative outer parts of the AD, which 
dominate the optical emission. The above effects are 
likely insignificant in this outer part of the AD. Thus, 
the redistribution of the AD radiation by these effects 
will not affect the measurement of 77, as long as we mea- 
sure the total SED, irrespective of its exact production 
mechanism. 

Here we derive M based on relatively sophisticated 
AD models, which include relativistic effects on the disk 
structure and photon propagation to the observer, and 
solve simultaneously for the vertical structure and radia- 



tive transfer of the disk. We apply the method to the PG 
quasar sample (jSchmidt fc Greenill983D. where Lboi is es 



timated bas ed on high quality opt ical (iNeugebauer et al 
1987[) . UV (IBaskin fc La;Mll2005l). far UV (IScotLeTa! 



20041 ). and soft X-rav (jBrandt et all 12000^ observations, 
and M is derive d based on high quality s pectroscopy of 
the H/3 region bv lBoroson fc GreenI ()1992D . The paper is 
organized as follows, in §2 we review the simple analytic 
derivation of M, and demonstrate that the AD L^, in 
the optical regime is rather well determined by the local 
blackbody AD models, and is only slightly modified by 
taking into account the vertical disk structure. We also 
show that the optical L^, is only weakly dependent on 
the radial disk structure, as set by a*. We then derive 
M for our sample. In §3 we estimate Lboi, and combined 
with M use it to compute rj. We discuss the correlation, 
or lack thereof, of 77 with parameters of interest, particu- 
larly M . In §4 we discuss various systematic effects which 
can affect the value of 77 and the observed correlation, in 
particular the uncertainty in M, disk inclination, opti- 
cal thickness of the AD emission, self-illumination, fore- 
ground extinction, and mass outflows. We summarize 
our conclusions in §5. 

2. ESTIMATING THE ACCRETION RATE 
2.1. Mass Estimates 

The spectral based methods for computing M outlined 
in t j2.2l and l2.3l reauire M estimates. We consider two sets 
of estimates, Mblr and Mg-, based on the broad emission 
line widths and the M — a.^, correlation, respectively. 

The first method requires characteristic velocities and 
radii for the Broad Line Region (BLR ) . We use the lumi- 
nosity radius relation of iKaspi et al.l ()2005[) to compute 
the BLR radius i?BLR, but with Lopt = vL^, measured at 
4861 A inst ead of 5 100 A . Inser ting this relation into equa- 
ti on (5) ofj Kaspi et all (l2000f) and using the H/3 FWHM 
V (|Borosori fc Grecnl ll992n . we can compute 



M, 



BLR = 1.5 X WMqLI-^I^^U 



3000' 



(1) 



where Lopt,45 is Lopt/10^ erg s~\ and v^om = H/3 
FWHM/(3000 km s-i). The quantities Lopt,45, "3000, 
and Mblr are reported in Table 1. 

A second method relies on the tight correlation 
between M and the st e llar velocity dispersion cr* 
(jFerrarese fc MerritH [2000t IGebhardt""etal. 2000). The 
major difficulty with this application is that the quasar 
light dwarfs the emission from the rest of the galaxy, 
making it particularly challenging to measure a.^, . There- 
fore, (T:, is only available for a handful of sources, and the 
measurements are likely less robust than those of inactive 
gala xies. For 13 of our sources , we also use estima tes 
from IDasvra et all ^2007,) and IWolf fc SheinT^ (|200l to 
compute Mc with the iTremaine et al.l (j2002 ) relation, 
and these values are reported in the second column of 
Table 2. 

2.2. Analytic Method 

The flux per unit area emitted by a thin AD is 
(jShakura fc Sunvaev|[l973l ) 



F 



3 GMM 
8^ i?3 



fc{r,a) 



(2) 
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where R is the radius, and fc{r,a) is a diniension- 
less factor, typically of order unity, which takes into 
account both the no-torque inner boundary condition , 
and the relatiyistic effects dNovikov & T hornd 119731 : 
IPaee fc Thornd 119741: iRTffert fc HeroldlllQQSri . This fac- 
tor, which approaches unity when r ^ 1, depends on 
a* and on the dimensionless radius r = R/Rg where 
Rg = GM/c^. 

We assume the disk emits locally as a blackbody at 
the effective temperature Tcff = (F/cr)^/'*. The full disk 
spectrum can now be obtained by integrating over the 
disk surface 



Btt^ cos iRg 



(3) 



where L^, is the observed specific intensity assumed to 
be emitted over 4tt steradians, i is the inclination to the 
line-of-sight, and is the Planck function. 

Since we are also primarily interested in the emission 
at relatively large radius, we will ignore the /c(r, a) de- 
pendence and use Toff = Tor"^/** with 



Tn = 



3c^ 



1/4 



Afi/2 ■ 



(4) 



Defining x = hu / [k^Tcti) and switching integration vari- 
ables yields 



where 



160 /GTr^/iG^c^ 



'^{v, Tin, To 



1/3 



cosi, (5) 



dx 



^8/3-1 



exp(a;) — 1 



(6) 



In the following, we will assume a constant Q w 1.93. 
This is approximately correct when j-'in ^ ^ t'out, 
where Vi-a_,ont are the frequencies of the peak of the emis- 
sion at Tin, out- For most of the PG quasars, this is a 
reasonable approximation as long as we evaluate L^^ in 
the visible range, since I'out is expected to be in the IR, 
while Via typically in the EUV. 

Defining Lopt = ^^-^i^ at 4861 A and solving equation 
© for M, we find 



M = 2.6Mo yr^i 



iopt,45 

cos i 



3/2 



(7) 



where Afs = A//10^ and iopt,45 = Lopt/lO^^ erg s~i 
Thus, the absolute accretion rate, AT, can be estimated 
from the observed Lopt, if M and i are known. 

2.3. Modeling Method 

In the previous section, we computed analytic expres- 
sions for Af from a simple, local blackbody model of 
an AD using several approximations. We now outline 
a calculation which takes advantage of more sop histi- 
cated spectra l models. We use KERRTRANS (|Agoll 
Il997t see also lDexter fc Agol|[2009l ) to calculate the disk 
integrated spectrum from a fully relativistic disk model 
(|Novikov fc Thornd Il973[ ). In this model the effects of 
varying the BH spin are included, with the assumption 
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Fig. 1. — Comparison of SEDs for local blackbody AD models, 
with different values of a* and M. All models have M = IMq yr~^ 
and cos j = 0.8. Note that the disk becomes hotter with increasing 
a» , but -Lopt (at the vertical dashed line) is generally independent 
of a, (excluding the a* = Af = 10^ model, which is too cold to 
be a viable AGN SED anyhow). Higher M AD models are colder, 
for a fixed M, as the SED i^pcak oc M~^/'^ , and are optically more 
luminous, as Lopt oc A/^/"^. Thus, if M is known, then M is set by 
Lopt ■ 

made that the innermost radius of the disk corresponds 
to the radius of marginal stability for circular orbits r,ns , 
and that no torque is present at this radius. We further 
assume no emission from inside Tms- 

We consider models in which the local emission at the 
disk surface is calculated with two different methods. 
It is either assumed to be a blackbody or it is com- 
puted using stellar- atmosphere-like calculations of the 
disk vertical structure. The latter (hereafter referred 
to as TLUSTY mod els) are essent i ally eq uivalent to the 
models described in lHubeny et al .' ('2000' fl and ar e com- 
puted using the TLUSTY code filubcnv fc Lanzj [T995). 

The relativistic model has four parameters: Af, 
a*, M and i. The spectra based on the full verti- 
cal structure calculations also require a choice of ass 
(jShakura fc Sunvaevlll973[ ) to determine the disk surface 
density. We assume ass = 0.01 for all models unless 
otherwise specified. 

These parameters completely determine the emission 
at all frequencies, and can be compared directly with the 
observations. Specifically, we will focus on matching to 
the observed luminosity at 486lA. We use the mass es- 
timates discussed in Wl.W to specify Af. We still need 
to specify i and a^,. For all but the highest masses, the 
inclination has the strongest effect on the derived ac- 
cretion rate, primarily due to the cosz dependence of 
the projected disk area. It is reasonable to assume that 
cos? 0.5 — 1 as nearly edge on systems are likely ob- 
scured. Thus, we adopt cosz = 0.8 as a characteristic 
value. We examine the implications of inclination de- 
pendence in later sections. 

Figure [1] shows the effect that varying a* has on rela- 
tivistic blackbody models. At fixed M and Af, the peak 

* The physical content of the models is identical. The only 
difference is that the spectrum for each annulus is computed via 
interpolation on a precomputed table of annuli, as described in 
Davis fc Hube nv (2006). 
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Fig. 2. — Comparison of local blackbody AD SEDs with the 
detailed atmospheric model TLUSTY, for different values of a* 
and a. All models have M = O.IMq yv~^, M = IO^Mq, and 
cosi = 0.8. Note the increasing atmospheric effects in the hotter 
parts of the disk, which are sensitive to the unknown viscosity 
mechanism (although the effect of a appears small). However, 
Lopt is rather insensitive to the atmospheric structure and radiative 
transfer effects, as it originates from colder parts of the AD, which 
is expected to emit locally close to a blackbody. Thus, Lopt should 
be mostly set by MM (see eq. [Sjl. 

of the spectrum increases with a*. Since M is fixed, 
Lboi increases with a* due to the increase in efficiency. 
However, this increase in Lboi mostly manifests itself as 
increased emission at high frequencies. At visible fre- 
quencies, the change with spin is much more modest. At 
IQ^Mq varying a* has a negligible effect. For IO^Mq 
there is a more pronounced effect because the optical 
photons are typically emitted much closer to r^s than 
at lower M, but the variation in predicted iopt from 
a, = - 0.998 is only ~ 50%. 

Figure [2] shows the differences between TLUSTY mod- 
els (ass = 0.01 — 1.0), and blackbody models computed 
with the same M and M for a, = and 0.9. The 
TLUSTY calculations tend to produce spectra which are 
nearly blackbody in the optical, but there is still a con- 
tribution to Lopt from the low energy tails of somewhat 
hotter annuli. These tend to be modified blackbodies 
(but with imprints of the Balmer edge) due to electron 
scattering, which shifts power to higher frequencies, leav- 
ing less flux at optical frequencies relative to the black- 
body prediction. At UV frequencies, this effect is more 
pronounced for higher a* , but the decrement is relatively 
independent of spin in the optical emission. 

For this M and M, the choice of ass has very little im- 
pact at any frequency, although the differences are larger 
at higher frequencies. This is due to the fact that the 
hot test, inner an nuli tend to be the lowest surface density 
in alShakura fc Sunvacvl (|1973f ) disk model. As discussed 
in lDone fc Davis (,2008i) , the spectra are generally insen- 
sitive to surface density as long as it is sufficiently large 
so that the disk is very optically thick. For higher Ed- 
dington ratio models, the surface density is low enough 
that ass can have a significant impact on the UV spec- 
trum (|Hubenv et al.,.20G 0). However, for the Eddington 
ranges of interest, the optically emitting annuli in AGN 
disks have large enough surface density that ass has very 
little impact on the optical spectrum. 



Fig. 3. — The derived M (in Mq yr -'^) as a function of M for the 
80 PG QSO in our sample. The a* = 0.9 and cos j = 0.8 TLUSTY 
model is used to derive M. The open and filled squares are derived 
using A/blr a-nd Ma-, respectively. The typical discrepancy in M, 
and in the implied AI, is a factor of 2 — 3, but can reach a factor 
of 10. Interestingly, the derived Af is only weakly correlated with 
M, although the AGN luminosity rises more steeply with M in our 
sample (see below). 

2.4. Accretion Rates 

With a*, i, M, and ass specified, estimating M is 
straightforward. We first calculate Lopt for models with 
the above parameters and different M. We then linearly 
interpolate to find the M for which the model optical flux 
would match the observed value. The resulting rates for 
our sample of PG quasars are shown in Figure |3] and 
reported in Table 1 and Table 2. A TLUSTY model 
with a^, — 0.9 and cosi — 0.8 was used to obtain these 
estimates. 

The M estimates are sensitive to the assumed M and 
the differences between Mblr and Ala- can be significant. 
Consistent with equation (O, a larger M, yields a smaller 
M and vice versa. The factors of 2-3 discrepancies in M 
yield comparable uncertainties in M. 

Figure m shows how M depends on a* and the method 
used to compute the surface emission. For each QSO in 
our sample, we compute five values of M using different 
models for the SEDs. The base model corresponds to the 
a* = 0.9 TLUSTY model used to derive the accretion 
rates plotted in Figure [31 The other four M estimates 
are used to compute ratios with the base model M in 
the denominator. The TLUSTY model with a^, = 0, 
generally provide a higher M. This is a 10% effect at 
low M, but can reach 40% at the highest AI. Blackbody 
models with a* = also tend to give similarly higher M 
for larger M, but can be 20% lower at low AI. Blackbody 
models with higher spin, a* = 0.9 and a* = 0.998, give 
about 10-20% lower AI for almost all AI. 

The objects with the two lowest masses correspond 
to very high Eddington ratios, which creates difficul- 
ties for the TLUSTY models. This is due to the fact 
that the annuli computed directly with TLUSTY do not 
cover the parameter range needed for the spectral mod- 
els, and extrapolation (rather than interpolation) is used 

to construct the spectra. As a result, the AI obtained 
with TLUSTY is probably an underestimate, resulting 
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Fig. 4. — The effect of different assumptions about a, and the 
local disk spectrum on the derived M. For each source in our 
sample, we compute M as described in i|2.3l using five different disk 
models. All values of M are plotted relative to the a* = 0.9, the 
TLUSTY based model results shown in Figure |3] The uncertainty 
in M is generally well below 40%, and is negligible compared to 
the errors resulting from the uncertainty in M. 

in the higher ratio of the blackbody models relative to 
this model that can seen in Figure HI 

The differences in M derived using different a* are con- 
sistent with Figure [1] in that higher a, yield larger Lopt 
when M is large, but comparable iopt when M is small. 
A larger model Lopt means that the observed Lopt can 

be matched with a lower M. The difference between 
TLUSTY and blackbody models are similarly consistent 
with Figure [2] in that TLUSTY models generally give 
lower Lopt, and thus require a higher M to match obser- 
vations. The uncertainty in a*, coupled with differences 
in the spectral models, corresponds to an overall uncer- 
tainty in M of ~ 20% at M = 3 X IO^Mq and - 40% 
at M = 3 X IO^Mq. We will find that this is signifi- 
cantly lower than the uncertainty associated with the M 
estimates. 

We are now in a position to evaluate how well our 
simple analytic relation ( §2.2|) approximates our more 
sophisticated fitting method ( §2.3|) . Treating the log- 
arithms of either (Lopt, M) or (Lopt, V3000) as sets of 
independent variables, we perform a linear least-squares 
fit to our M estimates. Our best fit relations are 

M = 3.5 Mq yr-i M^'-""' LI'^,^,,, (8) 

and 

M = 2.5 Mo yr-i v^oVo^'L^Z,, (9) 

Equation [8] is useful when Ma- is available, and equation 
[9] when only Mblr is available. 

As shown in figure O equation [9] provides a very good 
fit to our model-based M estimates, with an RMS error of 
AlogM — 0.024. Such precise agreement suggests that 
equation ([8]) or (jH) could be used in place of a detailed 
model fitting method for future work. A comparison 
of equations ([7]) and ([8]) demonstrates reasonable agree- 
ment between the simple analytic model and the best-fit 
relations, although the dependence on M is somewhat 
flatter than naively expected. This slight discrepancy is 
due to the effects of the inner boundary, which becomes 
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Fig. 5. — The model fit derived M plotted versus the best fit ana- 
lytic relations {eq[9]l using Lopt and i^aooo as independent variables. 
The RMS deviation of M from the analytic fit is AlogM = 0.024. 
This analytic relation provides a simple and accurate estimate for 
M, which does not require model fitting. 

more pronounced as M increases and the radius of opti- 
cal emission gets closer to the inner boundary. 

3. ESTIMATING THE RADIATIVE EFFICIENCY 
3.1. Bolometric Luminosity Estimates 

Estimation of the radiative efficiency of an accre- 
tion flow clearly requires a reliable measurement of the 
bolometric luminosity which is radiated. Since broad- 
band constraints on the SED are a priority, we focus 
on a sample of 80 relatively well observed P G quasars 
(|Boroson fc Greenlll993lBaskin fc Laoiil20050 . for which 
optical, UV, and X-ray data are available. However, the 
SED in the extreme UV, where the AD spectrum is ex- 
pected to peak, remain unknown. 

We compile data on the SED from several sources. 
The redshift, fluxes at 1549 A and 4861 A, the power- 
law slopes between these w avelengths a„uv, and t he H/3 
FWHM are aU taken from iBaskin fc Laori (|2005[ ). We 
also use a UV spectral slope auv computed between 
2000 A and 1400 A by A. Baskin (private communi- 
cation). Since these values are not published elsewhere, 
we report them in the fourth column of Table 1. If avail- 
able, the far UV slopes ctfuv and 1000 A flux es are taken 
from analysis of FUSE data presented in ' Scott et al.l 
(2004) and iShang eFall ((20050 The soft X-ray slope 
Qfx (0.2-2 keV) is computed from the H/? FWHM using 
the relation i n iBran dt et al.. (,2000.) . The flux at 1 keV 
is take n from iLaor ^ fc Bcharj (I20p8f). who tabula t ed th e 
data of IBrandt etTll ((2000) and lLaor fc Brandtl (|2002l ). 
All power-law slopes described in this text are (i.e. 
Fi, oc ly"). Fluxes are con verted to luminosities using 
the redshifts listed in .Baskin fc Laoii (|2005f ) and assum- 
ing a ACDM cosmology with Hq = 70 kms~^ Mpc"^, 
rim — 0.3, and = 0.7. 

^ F or sources lis ted in bo th Scott et al.l 1)20041 ) and IShang et al.l 
ifSOOSD . we use the IShang et a l. ( 2005) values. 
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Fig. 6. — Schematic piecewise power law SED used to estimate Lbol in our sample. In the loft panel we plot representative case A 
(dashed), case B (dotted), and case C (solid) model SEDs as described in the text. The dark shaded areas denote frequency ranges where 
the SEE) is computed using observations for individual models, while the light shaded areas show frequencies where we have assumed 
characteristic slopes motivated by QSO observations at these frequencies. The unshaded regions are interpolations or extrapolations. Note 
that the optical to UV frequency range is expanded relative to X-ray for clarity. The right panel is a 'zoomed-in' plot of the UV to soft 
X-ray region of the left panel to highlight differences between the three cases of model SEDs. The three different possible cases for the 
SED interpolations imply a typical (RMS) uncertainty of 0.24 in logLjjol- 



With these data we compute Lboi using a piece-wise 
power law representation for the SED. To clarify the 
discussion that follows, we plot a characteristic exam- 
ple SED in Figure [5] It is relatively straightforward 
to estimate the continuum for wavelengths longer than 
1549 A and shorter than 62 A (0.2 keV). First, we 
exclude the observed infrared bump shortward of 1 /im. 
This emission is thought to be reradiated by dust and we 
are only interested in the direct emission from the AD. 
Therefore, we assume power-laws with a = 1/3 below 1 
fim and a = —0.3 between 1 /im and 4861 A. Between 
4861 A and 1549 A, we assume a power law slope equal 
to the measured aouv The spectrum between 0.2 and 2 
keV is assumed to have a power-law slope equal to 
and is normalized to match the 1 keV flux. From 2 keV 
to 50 keV, we assume a = —1, and above 50 keV we 
assume a — —1.5. 

We consider three different prescriptions for the unob- 
served extreme UV spectrum. In case A, we assume a 
power-law with a = —1 between 1549 A and 1000 A, 
and another power-law is fit through 1000 A and 0.2 
keV. In case B, we use aouv for the power law slope be- 
tween 4861 A and 2000 A, but use the measured a^v 
from 2000 A to 1000 A. Between 1000 A and 200 A, we 
assume a = ctuv ~ 0.307, and a power-law is fit between 
200 A and 0.2 keV. Finally, for those sources with FUSE 
data, we also consider a case C. We fit a power law be- 
tween the 1549 A and 1000 A fluxes. From 1000 A to 
200 A, we use afuv, and a power-law is again fit between 
200 A and 0.2 keV. Note that Aa = -0.307 in case B 
corresponds to the mean difference between a^v and afuv 
in those sources where afuv is available from FUSE data. 

The resulting SEDs for case A, and for B or C are plot- 
ted for each source in Figure [T] The plots are arranged 
in order of increasing M. We also show an a* = 0.9 
relativistic blackbody model with the known A/, and a 
value of M which has been fit to match the luminosity 
at 4861 A ( see ^2.3^ . In most of the quasars, the SED 
peaks in the observable far UV (rather than the unob- 



servable extreme UV). The X-rays generally contribute 
only a modest fraction of the overall luminosity. 

For the majority of sources there is reasonable agree- 
ment between the various SED models, but there are 
several sources with rather discrepant results. For cases 
B and C, there are a number of sources where a^v or 
afuv seem to be rising or falling much more steeply than 
appears likely. This is presumably due to some combi- 
nation of slope and fiux measurement errors, dust red- 
dening, and variability (the X-ray and UV fiuxes are not 
contemporaneous). In sources with FUSE observations 
there is sometimes a mismatch between the luminosity 
estimates at 1000 A and the lBaskin fc Laoil ()2005[ ) data 
at 1549 A, presumably due to instrumental uncertainties 
and variability. 

The case A SED, constructed to form a smooth tran- 
sition between the UV and soft X-rays, often seems to 
be the "most reasonable" . However, it is important to 
avoid enforcing a universal SED preconception, and so 
we adopt the approach of using as much UV data as 
possible, even when it gives apparently unlikely FUV ex- 
trapolations. Therefore, in sources where FUSE data is 
available, we use the case C SED to estimate Lboi- In 
sources where it is not available, we use case B. We re- 
port the resulting Lboi estimates in the fifth column of 
Table 1. We find that this choice increases the scatter in 
our ?7 estimates over what we would have obtained us- 
ing case A, but does not significantly change the overall 
trend with M. We make one exception for PG 0049-1-171, 
which has afuv = 4.1. Such an unphysically steep slope 
leads to an extreme overestimate of Lboi, so we neglect 
it in all further analysis. 

There are several potential sources of uncertainty in 
the Lboi estimate derived from these SEDs, but in the 
vast majority of sources the dominant uncertainty is the 
far-UV extrapolation. To estimate this uncertainty, we 
assume Aa = ±1 for far-UV extrapolation afuv (which 
is equal to auv — 0.307 for case B) and compute the re- 
sulting range of £boi- This is considerably greater than 
the typical measurement errors in the sources where afuv 
is measured directly with FUSE. We adopt a larger un- 
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Fig. 7. — The SED derived for each object. The sohd curve is the SED used to compute Lbol- It corresponds to case C for sources with 
FUSE data and case B otherwise. The dashed curve shows the case A SED for comparison. These SEDs are described in detail in ^3.11 
The dotted curve is the a* = 0.9 local blackbody AD model with M = Mblr that matches the SEDs at 4861 A. The objects are ordered 
by ascending order of M. Note that the model fit is hotter than the observed SED for low M objects, and colder for the high M objects. 
Thus, the observed AGN SED in our sample is inconsistent with a fixed a, AD model 




Fig. [7] (continued) 
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certainty to conservatively account for any systematic 
errors, possible effects of dust reddening, and potential 
complexity in the unobserved part of the SED. Since this 
generally yields a large uncertainty that dominates other 
sources of error (e.g optical or X-ray variability) we as- 
sume the contribution from other sources is negligible. 

Note that in Figure[71 the AD model, which has a fixed 
77 = 0.16 (a* — 0.9), systematically overpredicts the ob- 
served FUV SED in the lowest M objects and produces 
a higher overall Lboi- In contrast, the AD model sys- 
tematically underpredicts the observed FUV SED for the 
highest M objects, and produces a lower total £boi- This 
systematic trend suggests that a single value of r\ will not 
be consistent with all objects, but instead implies that 
77 needs to increase from low to high M, as we discuss 
below. 

3.2. Radiative Efficiencies 

Now that we have estimates for Lboi and M we can 
estimate r/ using 

Lho\ 



V 



(10) 



The factor of cos? accounts for the inclination depen- 
dence of the observed Lbob if it originates in a thin AD0 
Here we assume cosi = 0.8 as noted above (§2.3). 

A linear least squares fit for log 77 as a function of the 
logarithms of M, Lopt, assuming a constant error for each 
77 estimate, yields best fits of 



77 - 0M3LU%L-X45K 



fO.89 



or 



77 = 0.086iU%L-°«6..i.78 



Dl,46-^opt,45''^3000' 



(11) 



(12) 



with rms deviations of Alogr/ = 0.024 from the true 
values. A comparison of these best-fit values with our 
TLUSTY based model estimates is shown in Figure [HI 
We would have arrived at nearly identical relations if we 
had simply inserted equations ([8|) and ([9]) into equation 
pUj) . The best fitting relation when Lboi is not available 

77 = 0.08L,-p",l|M°■«^ (13) 

with an rms deviation of A log 77 = 0.26. This is still 
interesting, but considerably poorer than the three pa- 
rameter fit when Lboi is available. 

To understand these best-fit relations we look at the 
expected relation for the simple Newtonian AD case. In- 
serting equation ([7]) into (fT0|) yields 

77 = 0.068Lbol,46iopM5^8(cOSi)'/', (14) 

or using equation ([T]) 

r, = 0.11Lboi,46iopM5^3'ooo(cos i)'/'- (15) 

The assumption that cosi ^ 0.5 — 1, combined with the 
(cosi)^/^ dependence suggest that inclination uncertain- 
ties are not likely to be a significant source of uncertainty 

® For a flat Newtonian AD there is an additional factor of 2 in 
the numerator. However, relativistic beaming and electron scat- 
tering induced limb darkening modify the angular distribution at 
higher frequencies, invalidating the Newtonian approximation. We 
therefore neglect the factor 2 for simplicity. 
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Fig. 8. — T he derived rj plotted versus the best fit analytic rela- 
tions (eq. II2I I using Lopt, Lboli a-nd '^aooo a-s independent variables. 
The RMS deviation of r] from the analytic fit is A log 7] = 0.024, 
i.e. an accuracy of 5%. 

in rj. With the inclination dependence removed, these 
results are in approximate agreement with our best fit 
relations, confirming that simple Newtonian blackbody 
model captures most of the relevant physical effects. 
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Fig. 9. — A plot of r] versus M for the 80 PG quasars in our 
sample. The symbols correspond to estimates mad e with Mblr 
(blue, open squares) and Mp- estimates obtained bv lDasvra et al] 
teoOTi) (filled squares) and I Wolf fc SheinisI 12001) (filled triangles). 
The dashed solid lines connect the two sets of M — estimates for 
the 13 sources for which both A'/blh, a-nd A/o- are available. The 
rj error bars account for uncertainties in L^oi only. The horizontal 
dotted lines corre spond to 17 = 0.038, 0.057 and 0.31, the theoretical 
efficiencies from INovikov fc Thornel {1973) for a, = —1, 0, and 
0.998, respectively. The thick solid line in the lower right hand 
corner of the plot shows the displacement r) — M plane that would 
occur for a 0.4 dex shift in M. Two objects, PG 2209-f 184 and 
PG 1512-f307, have r; > 1, although their errors are consistent 
with 77 ~ 1. Note that most objects fall within the theoretically 
expected range for t). The outlying objects are within the expected 
uncertainty in M. There is a clear rise in r] with M , with a best 
fit relation r) = 0.089M0-52. 

In Fig ure m we plot 77 as a function of M using M blr, 
and M„ (IDas vra et al']l2007t IWolf fc Shei^ 120081) . We 
report the BLR and a based estimates of 77 in the last 
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columns of Table 1 and Table 2, respectively. For the 
13 objects with both M„ and Mblr estimates, ry is plot- 
ted for both cases and the symbols are connected with 
dashed lines. The rj error bars are computed from the 
uncertainties in Lboi- The uncertainty associated with 
the value of M is likely larger for most sources, but is 
not included. Since the uncertainty in M affects both 
the abscissa and ordinate it will shift points diagonally 
in the plot (as seen for the dashed lines). The thick solid 
line in the lower right hand corner of the plot shows the 
displacement in the rj — M plane that would occur for 
a 0.4 dex shift in M . The impact of this uncertainty is 
discussed in detail in ij4.2l 

The mean values of rj are log?] = —1.05 ± 0.52 and 
—0.91 ±0.30 for the estimates made with A/blr and M^. 
These values are consistent with the values derived from 
the Soltan argument (see §1) which is based on com- 
pletely independent arguments. In addition, we find a 
clear correlation of rj and M . We find a Spearman rank 
correlation coefficient r = 0.85, with a significance (prob- 
ability) Pr — 1 y. 10^^^ and a best fit relation 



Tj = 0.089Af; 



0.52 



(16) 



The majority of the sources are consisten t with 0.057 < 
t; < 0.321, the iNovikov fc Thornd ([1973') efficiencies for 
< a* < 0.998. At high M, several sources have -q > 
0.321, but are generally also consistent with ij < 0.321. 
There are two sources with r] > I (PG 2209-1-184 at 
logM = 8.2 and PG 1512-^307 at log A/ = 9.2) and Fig- 
ure [7] shows that both have two of the most steeply rising 
FUV slopes in the sample. These FUV slopes are diffi- 
cult to reconcile with the X-ray flux level and are likely 
erroneous, leading to overestimation of iboi and, there- 
fore, overestimation of r/. In contrast, At low Af , there 
are numerous sources with relatively small error bars, 
based on the Lboi estimates, for which rj is signiflcantly 
less than 0.057. Our result is consistent with the result 
of Collin et al. (2002), who assumed a constant 77 model, 
and derived an increasing Lboi/iEdd with decreasing Af, 
with iboi/iEdd reaching 1-100 for an assumed ij = 0.32. 
Assuming iboi/^Edd is at most only slightly greater than 
unity, as inferred for our sample, this implies 77 ^ 0.01 
for the low M objects in their sample, broadly consistent 
with our result. 

In Figures [To] and [TTl we plot the variation of the A/blr 
based rj with Lboi/^Edd and the radio loudnes s parame- 
ter R = /6cm//44oo from ILaor fc Beharl (|2008n . In each 
panel we provide r and Pr for each distribution. Fig- 
ures [10] shows a modest (r = 0.54) correlation of 77 with 
Lhoi/ L^dd- In Figure [TT] there is a tendency for radio 
loud quasars to have higher average 77 than radio quiet 
sources, although the distributions overlap and the re- 
sulting correlation is weak (r = 0.25). Given that Figure 
[9] shows a correlation of 77 with M, such a tendency is 
not unexpected since there is a well-known correlation of 
R with M (e.g. [Lao r 2000). This is in qualit ative agree- 
ment with the results of ,Bian fc Zhaol ()2003D of a higher 
77 for the radio loud quasars. 

3.3. iopt/fjboi Ratio and its Implications 

In Figure [T^ we plot the ratios Lopt /iboi versus the 
FWHM of the H/3 line. The observed ratio of Lopt/iboi 
shows almost no trend with the FWHM. The distribution 
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Fig. 10. — Correlation of rj with Eddington ratio. The Spearman 
rank correlation coefficient r, and corresponding significance Pr 
confirm that a modest correlation exists. However, this correlation 
is likely induced by the stronger correlation of rj and M {Fig|9]l, 
and the correlation of the Eddington ratio with M (see below). 
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Fig. 11. — Correlation of t] with radio loudness, R. The Spear- 
man rank correlation coefficient r, and corresponding significance 
Pr indicate that only a weak correlation exists. This correlation 
is likely induced by the correlation of rj and M, and the known 
correlation of radio loudness and M in the PG quasars. 

of log iopt/ibol has a mean of -0.96 with a standard devi- 
ation of 0.27. The errors are generally larger for objects 
with low iopt/iboi- In these objects, the far UV tends to 
account for a larger fraction of Lboh so the uncertainty 
associated with the extrapolation into the extreme UV 
and soft X-rays tends to be more important, giving a 
larger overall uncertainty in Lboi- When we instead use 
the more uniform case A SED (not plotted) , we find con- 
siderably lower scatter, but the absence of a trend with 
FWHM remains. 
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Fig. 12.— Ratio of Lopt/ibol as a function of the FWHM of H/3 
for the objects in our sample. The errors are computed assuming 
an uncertainty of Ao = ±1 for ctfuy for all objects. The dashed 
line is a simple analytic estimate (eq. I17| l for constant rj. The open 
squares are also constant -q estimates, but instead based directly 
on our M estimates, as described in i|3.3l The ratio of iopt/^bol 
shows a relatively small spread, and no correlation with the FWHM 
of H/3, and demonstrates again that constant r} thin AD models are 
excluded. 

The strong correlation of rj and M is consistent with 
the lack of a correlation between the ratio iopt/iboi and 
FWHM in Figure [l2l Solving equation (fTSj) for Lboi, we 
obtain 



L 



opt 



0.09877o:!l°p\%^^3'ooo, 



L 



bol 



(17) 



where 770.1 — 77/O.I and we have assumed cosi = 0.8. 
Therefore, if 77 were constant, this ratio should vary 
roughly as the square of the FWHM with only a very 
weak dependence on I/o^t- We plot this relation as a 
dashed curve in Figure [12] For simplicity, we ignore the 
Lopt dependence and set Lopt,45 — 1 in equation (IT7)) . 
Alternatively, we can use our M estimates derived from 
model fitting to speciiy Lboi via equation PH)) assuming 
constant rj = 0.089 and cosi = 0.8. These estimates 
are plotted as open squares in Figure [12] and follows the 
dashed curve with little scatter. This shows that the 
evolution of Lopt / Lboi derived from the TLUSTY mod- 
els is largely captured by our simple analytic estimate. 
In either case, the expected evolution of Lopt/Lboi for 
constant rj is not observed. 

The above relation assumes that M is simply a func- 
tion of Lopt and FWHM. Since there is a potential for 
systematic uncertainties in these M estimates, it is use- 
ful to consider what relationship between Lboi and M 
would be needed to reconcile a nearly constant Lopt / Lboi 
with a constant rj. From equation ([T^. we see that 

1 /2 

M oc Lj^^j or, equivalently, Lboi/LEdd cx M is required to 
keep 77 constant. The point like selection criterion for the 
PG sample selects against low Lboi/LEdd objects, where 
the host galaxy becomes dominant. This likely induces 
the relatively high values and small spread in Lboi/LEdd 
(Figure [TS)) . yielding Lboi oc M. Thus 77 oc M^-^, as ob- 
served, results from the small spread in both Lboi/LEdd 
and Lopt /Lboi- While the first effect is possibly a se- 
lection effect of the PG sample, the second is almost 
certainly not. 

We note that a relatively uniform ratio of Lopt /Lboi 




Fig. 13. — Ratio of Lboi/iEdd versus M for all 80 PG 
quasars in our sample. The symbols correspond to estimates 
m ade with Mblr (blue , open squares) and Mg estimates obtained 
by Dasvra et al.l IpOOTl ') (filled squares) and IWolf & Sheinisl 1(2003 ') 
(filled triangles). The y error bars account for uncertainties in Lboi 
only. The deficiency of objects with Lbol/^Edd < 0.1 likely results 
from the selection of point-like objects for the PG sample, as ^Edd 
is proportional to the host bulge luminosity. The rise of the lower 
limit to LbolZ-^^Edd in the lowest M objects may be due to M es- 
timate errors, if all low M objects have a higher M . Or, it could 
also be a selection effect, if low M objects reside in disk galaxies, 
where the AGN needs to outshine the disk light as well. 



has b een found by previo us authors (e.g. [Elvis et all 
fl99l [Richards et alll2006l ). and is, in fact, frequently 
assumed in the literature to estimate Lboi when only op- 
tical constraints are available. Using simple templates 
with fixed optical-to-UV SEP, b ut accounting for X-ray 
variation, [Marconi et al.[ ()2004[ ) find weak evolution in 
the bolometric correction that corresponds to an increase 
in Lopt /Lboi with M . However, this is still a much weaker 
dependence than expected for constant Lboi/LEdd and 77. 
Studies with large QSO samples that assume a constant 
bolometric correction and BLR mass est imates generally 
find Lhoi/LEdri either nearly constant ([KoUmeier et all 
[20061 ) or decreasing slightly jStei nhardt fc Elv id 120101) 
with M . As with our PG sample, reconciling these ob- 
servations with constant 77 not only requires the BLR 
estimates to be in error, but Lboi/LEdd would have to 
scale nearly linearly with M . 

4. DISCUSSION 

In the previous sections we outlined and implemented 
a prescription for estimating M and, ultimately, 77 by as- 
suming ADs are radiatively efficient. Our mean 77 iscon- 
sistent with con straints derived using the [Soltad ([1982[) 
argument (e.g . 'Yu & Tremainc 2002; El vis et al.[[2002[: 
[Marconi et all [2004: Bargcr et al. 2005), and with the 
theoretically expected range of rj for thin AD. This sug- 
gests that the radiatively efficient, thin disk is a viable 
model for accretion flows in QSOs, at least at the radii 
where Lopt is generated. However, what drives the ob- 
served strong correlation of 77 with Af ? We consider three 
basic possibilities: the rj — M correlation is real, possibly 
due to a correlation of a* with M; actual accretion flows 
differ from the radiative efficient model in an M depen- 
dent manner; or the input parameters {M or Lboi) are 
not estimated reliably. We discuss each of these possibil- 
ities in the following subsections. 

4.1. Effects of error in Lopt o.nd Lboi 
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Our estimates of Lboi are subject to two primary un- 
certainties: the unknown inclination and our inability to 
observe in the the extreme UV. The unified model of 
AGNs suggests that our sample, which are all t ype 1, 
are viewed not far from face on (| Antonuccil 1 1 993*) . Evi- 
dence suggests that the opening angle increases with lu- 
minosity, but is < 60° at the higher lumino s ities charac- 
terist ic of our sample (jPolletta et alj 120081 : iReves et alj 
|2008() . Therefore we expect our sample to cover a range 
of cosi ^ 0.5 — 1 (or smaller). To the extent that the 
intrinsic emission is isotropic, the geometric dependence 
is simply Loptj -^boi oc cosi. As noted in H'6.2[ this results 
in an (cosiY^'^ dependence for rj, which is only 0.15 dex 
of scatter for a factor of 2 uncertainty in cosi. However, 
this neglects effects due to atmospheric limb darkening 
and relativistic beaming, which are wavelength depen- 
dent and, therefore, give different i dependence for Lopt 
and Lboi- For the relevant range of M and Lboi, an ex- 
amination of the model SEDs suggests a scatter of < 0.3 
dex in rj at fixed M, for a uniform distribution in cosi 
between 0.5 and 1. 

Our estimates for Lboi and its uncertainty are discussed 
in The uncertainties in Lboi, which determine the 

plotted error bars for rj, entirely reflect the uncertainties 
in our FUV slope extrapolations. This uncertainty could 
account for much of the scatter in the rj — M correlation 
(Figure |9]), but not the overall trend. A simple examina- 
tion of Figure [7] shows why the rj — M correlation cannot 
be easily explained by "hidden" emission in the EUV. At 
the high M end the models predict a rollover in the spec- 
trum in the observable FUV, which is not seen in most 
sources. Thus, any unobserved emission in the EUV, 
only makes the discrepancy larger. At the low M end, 
the models continue rising into the EUV even though 
there is often a clear flattening or rollover in the SED 
in the observed FUV. Since the X-ray flux is typically 
below the FUV, one would require pathologically double 
peaked SEDs to provide the "hidden" emission in these 
cases. 

The fluxes used to construct our SEDs have been cor- 
rected for Galactic dust reddening and neutral absorp- 
tion, but our analysis does not account for intrinsic red- 
dening, which could result in a systematic underestimate 
of Lboi relative to Lopt and, therefore, an underestimate 
of 7]. An an alysis of a large Sloa n Digital Sky Survey 
QSO sample (|Hopkins et al.|[2004[ ) concluded reddening 
is modest in most sources, with only a small fraction 
which are reddened significantly (i.e. < 2% with E(B- 
V) > 0.1). However, some of the SEDs in Figure [7] give 
the impression that reddening might possibly be affecting 
the FUV and our extrapolation into the extreme ultravi- 
olet. Such sources tend to have a smaller fraction of their 
emission in the FUV and show up as the lowest 77 sources 
in their M bin. Since the error bars are determined by 
the uncertainty of the EUV to soft X-ray contribution to 
the SED, which is relatively small for these sources, they 
also tend to have relatively small error bars in Figure 
[9l These arguments suggest that some of the lowest ry 
estimates at low M are likely underestimates due to red- 
dening. However, it seems unlikely that reddening alone 
can explain the overall correlation of 77 with M, as this 
would require systematic reddening of all sources at low 
M, which is not evident in the SEDs. 



Contamination of optical emission from sources other 
than the AD could potentially lead to overestimates of 
M and underestimates of ij. Potential sources of such 
contamination are host galaxy or jet emission. However, 
significant optical jet emission is likely constrained to 
the small fraction of core dominated, blazar like radio 
loud quasars, which account to < 5% of the PG quasars. 
Furthermore, these tend to be at high M, where t] is 
already uncomfortably too high. The host galaxy must 
be present at some level and would have its greatest effect 
for low iboi and high M sources. I ndeed a significant 
host galaxy contribution is inferred bv lShen et al.l ()2010f ) 
for SDSS QSOs with Lopt < lO'*^ erg s'K If we adopt 
their correction for our Lopt, which is measured at 4861 
A (as opposed to 5100 A in their work) we find that 
several of our lowest r] sources at low M our shifted to 
higher values, although many still remain with r] < 0.038, 
the nominal lower value for maximally spinning, counter 
rotating BHs. 

4.2. Effects of error in M 

Our M estimates are described in H2.U The uncer- 
tainties on these values are somewhat difficult to esti- 
mate since they are subject to systematic errors which 
are not easily quantified. The Mrt. r estimates utilize 
the empirical i?BLR — Lopt relation of iKaspi et al.l (|2005[ ) 
which is calibr ated against a samp le of reverberation 
mapping radii (jPeterson et al.|[2004D with 40% intrinsic 
scatter. The reverberation mapped M estimates them- 
selves show_a^_scatter_of 2.6-2.9 relative the M — a rela- 
tion ()Onken et al.ll2004D . Therefore, a factor of ~ 3 (0.4 
dex is frequently stated) is widely reported as a charac- 
teristic uncertainty in these estimates. However, argu- 
ments have been made that suggest the true error could 
be either smaller or larger than this value. A thorough 
discussion of systematic errors, includi ng inclination de- 
pendence, is provided bv iKrolifl ()200lD . Lower mass er- 
rors (some ^0.2 dex) on BLR based estimates have also 
been suggested by some authors (Kollmcicr ct al. 200^ 
'Shen et al.l [200l ISteinhardt fcElvisI l2010t iKellv et all 
,2010), based primarily on the presences of sharp fea- 
tures in the inferred Eddington distributions. Also, the 
small transition range of log M = 8.5 — 9, where the PG 
quasars transform from radio- quiet to rad io-loud, sug- 
gests a small uncertainty in M ()LaoH[2000() . 

The M — a estimates are rather uncertain due to 
difficulties in measuring a in galaxies that host bright 
quasars. We note that the mean absolute deviation be- 
tween the two estimates is AlogAf — 0.6 for the 13 
sources with estimates (Fig. although it is not 
clear which is the more reliable estimator. 

Our model fitting estimates roughly give 77 oc M at a 
fixed luminosity (ea ll4l) . Therefore, errors in M and 77 
are correlated with AlogM ~ A log 7;, scattering data 
points diagonally in Figure [9l Since the inferred corre- 
lation is approximately rj oc M^/^, this scatter will have 
a significant projection onto the observed trend. If we 
make the reasonable assumptions that mass errors are 
approximately symmetric and centered on the correct 
M, we can consider the effect they will have on our r; 
estimates. 

In order to model the effects of mass uncertainties, we 
have implemented a simple Monte Carlo procedure to 
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Fig. 14. — Simulated distributions showing the effects of mass 
measurement errors on the derived -q. For comparison, our mea- 
sured rj — M distribution is overplotted as blue squares and the 
assumed mass-independent intrinsic r) is shown as a dashed hori- 
zontal line. The panels show distributions with different combina- 
tions of (Tint and (Torn the standard deviations on the log normal 
distributions modelling the intrinsic (true) mass distribution and 
the mass measurement error, respectively. To match the observed 
rj — M distribution, the PG sample needs to have both a narrow 
intrinsic M distribution, (Tint = 0.4, and a relatively large mea- 
surement errors, crcrr = 0.6. With decreasing (Torr and increasing 
(Tint the predicted r) — M correlation weakens (middle panel), and 
eventually disappears (lower panel). 

generate observed distributions of parameters based on 
our AD spectral models. The procedure assumes log nor- 
mal distributions for intrinsic mass Mint and Eddington 
ratio £ = iboi/-^Edd with means /Xjnt and standard 
deviations tiint and (T^, respectively. We further assume 
that the uncertainty in mass AM = Mobs — Mint follows a 
log normal distribution with a mean of zero and standard 
deviation CTcn-- Inclination i is assumed to be distributed 
uniformly in cosi between 0.5 and 1. We adopt a single, 
fixed value of r/int, as our focus is on determining whether 
the rj — M correlation could arise solely or primarily from 
errors in our mass estimates. 

Each Monte Carlo realization is generated by the fol- 
lowing steps: 1) We draw Mint, £, and «int randomly 
from the above distributions. 2) We compute Lboi,obs = 

^Mintcosi and Mint = £Mint/(?7intC^)- 3) We compute 
the observed Lopt from our AD spectral model using 
Mint, Mint, and iint- 4) We draw AM from a log nor- 
mal distribution and compute Mobs = Mint + AM. 5) 
We use the fitting procedure described in ij2.3l (with 
cosiobs = 0.8) to infer an observed accretion rate Afobs 
from Afobs and Lopt- 6) We infer an observed efficiency 

77obs using equation (|10l) with Afobs and iboi,obs- Note 
that steps 4 and 5 would simply invert steps 2 and 3 if 

Afint = Afobs and iint = iobs- 

This algorithm is repeated a large number of times 
to map out the distribution of Afobs and ?7obs for each 
set of input parameters /i^, /iAf, ai, Uint, CTcrr, and 
?/int- We focus on distributions with /iint = logM = 8, 
m =\og£= -0.4, ai = 0.4 dex, and ?7int = 0.08. The ob- 



served distribution of r] and M, along with our assump- 
tion that AAf is log normal with zero mean, strongly 
constrain these choices for /iint and ?7int. We choose /if 
and (T( to approximately reproduce the observed distri- 
bution iboi/^Edd seen in figure [T3l although the optimal 
values depend slightly on (Tint and CTqit- 

The resulting distributions are plotted as contours in 
figure I14[ with the observed rj — AI distribution overplot- 
ted for comparison. Each panel shows a different combi- 
nation of (Tint and (Terr- From top to bottom in figure fT4l 
we increase (Tint , while keeping (Tint + Ccn- — 1 • This con- 
straint is motivated by the observed range of M, which 
is roughly consistent with (Tint + <7err ~ 1 • 

We find that we can largely reproduce the observed 
r/ — M correlation by choosing a sufficiently small (Tint 
and a sufficiently large CTorr, although difficulties remain 
for this interpretation. As noted above, (Torr 0.4 dex is 
often asserted, but this value would be insufficient to re- 
produce all of the correlation by mass errors alone. Even 
if one allows that (Terr = 0.6 dex and CTint = 0.4 dex are 
theoretically viable, this distribution has mixed results 
simultaneously reproducing all of the observables. 
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Fig. 15. — Simulated distributions showing the effects of mass 
measurement errors on the observed distributions of ibol/^'Eddi 
Lopt, and Lopt/ibol- For comparison, observables from our PG 
sample are overplotted as blue squares. The distributions are com- 
puted assuming (Torr = 0.6 dex and (Tint = 0.4 dex, correspon ding 
to the rj — M distribution plotted in the top panel of figure 1141 
Note that although the observed distributions in the upper and 
lower panels can be reproduced, the observed I/opt vs. M corre- 
lation is stronger than expected if the error in M is as large as 
required. 

This can be seen in figure [121 which shows the two 
dimensional distributions of Mobs with fvboi/.^Edd, ^opt, 
and Lopt/f^boi for the (Torr = 0.6 dex and CTint = 0.4 dex 
distribution. Even though our intrinsic distribution of 
iboi/^Edd is mass independent, an anti-correlation with 
M arises which is qualitatively consistent with the ob- 
served distribution. Reasonable agreement is also found 
for fjopt/^boi- However a significant discrepancy remains 
between the simulated and observed Lopt distributions, 
which shows a much tighter correlation of Lopt with M 
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in the observed data. The breadth of the simulated dis- 
tribution is a result of the scatter induced by mass errors. 
We can get a narrower, more correlated distribution by 
reducing (Terr to 0.2 dex, but this produces a steeper de- 
pendence of Lopt on M and doesn't produces a large 
enough rj — M correlation. 

It is plausible that some of remaining discrepancy could 
be addressed with more complex distributions for the 
Mint, and AM. We have no strong empirical or the- 
oretical motivation to believe that the real distributions 
of the above variables are log normal. In fact, these as- 
sumptions lead to an A/obs distribution which is also log 
normal. This conflicts slightly with our observed M dis- 
tribution, which is somewhat "flatter" (i.e. has negative 
kurtosis) than a log normal distribution with the same 
variance. We have little empirical knowledge of the true 
distribution, which is subject to the non-trivial selection 
effects of the PG sample. Therefore, we can only hope 
to approximately quantify the effects of M uncertainties 
and we adopt the above assumptions largely for the sake 
of simpl icity. A mo r e sop histicated analysis, along the 
lines of IKellv et al.l (|2010( ) is beyond the scope of the 
present work and possibly not amenable to our limited 
sample size. 

Regardless of the precise details, this analysis relies 
on the (plausible) assumption of an approximately sym- 
metric mass measurement error and a centrally peaked 
intrinsic mass distribution. At high mass, such a dis- 
tribution could arise because large M objects are rare, 
but the turnover at low M would require selection effects 
that tend to exclude low M objects. This could be the 
case if low M systems are Eddington limited and tend to 
have too large of a host fraction to make it into our PG 
sample. 

We note that our treatment of inclination has little ef- 
fect on the resulting distributions. We model a uniform 
distribution of cosi from 0.5-1, but assume a single value 
of cos i — 0.8 for our observed inclination. In the absence 
of mass errors, this introduces a modest broadening of 
77obs about r^int , and the magnitude of the effect is consis- 
tent with our above estimates in M.ll However, this is 
dwarfed by the broadening introduced by any plausible 
mass error (i.e. any tJcrr > 0.2 dex). A dependence on 
inclination may also be present in our virial M estimates, 
since the translation of the line width to a virial veloc- 
ity will in general depend on the inclination distribution 
of the emitting gas, which may (in turn) be related to 
the inclination of the AD. Therefore, inclination uncer- 
tainties may have a much stronger impact through their 
effect on M than their effect on either Lopt or Lboi- How- 
ever, we consider this effect implicitly modelled by our 
mass error analysis above. 

In addition to the distributions shown in figures fT4l and 
1151 which do not model dust reddening, we have gener- 
ated Monte Carlo simulations that incorporate the effects 
of reddening on the observed Lboi and iopt. We consider 
several different redd ening curves, including SMC-like 
(jRichards et al."2003') and those of iCzernv et all (|2004f ) 
and Gaskell ct al. (2004). As expected, we find that red- 
dening tends to lower the ratio of iboi/-^opt, decreasing 
the inferred rj. This result is independent of reddening 
curve, although the amount of extinction required de- 
pends on the curve used. We find moderate amounts 
of extinction (E(B - V) < 0.1 for SMC-like reddening) 



are all that is needed to explain some of low rj sources 
that populated the lower envelope of the observed dis- 
tribution. The decrease in rj tends to be slightly larger 
for low A/, which have a larger fraction (relative to high 
M systems) of their bolometric output in the UV in our 
model SEDs. This effect alone introduces a slight posi- 
tive correlation between rj and M, but the magnitude of 
slope change is small and cannot contribute significantly 
to the observed trend for reasonable extinction values. 

Our conclusions from this analysis is that mass mea- 
surement errors can introduce a spurious correlation of 77 
with M, and that this effect almost certainly contributes 
to some of the inferred correlation. However, given the 
tight correlation between Lopt and M, and the require- 
ment of a very narrow mass distribution (crint ^ 0.4 dex), 
we disfavor M errors alone as an explanation for the ob- 
served correlation. 

4.3. The accretion disk model 

Since our rj estimates are all ba sed on the bare, ra- 
diativ ely efficient thin disk in odel (jShakura fc SunvaevI 
1973; Novikov fc Thornelll973l ). it is conceivable that the 
evolution in 77 with M is a result of inapplicability of 
the underlying spectral model. Indeed, it is well known 
that this model has difficult i es rep roducing AGN obser- 
vations. iKoratkar fc Blaej (|1999l) offer a thorough re- 
view of discrepancies between models and observations. 
In particular, it fails to account for the X-ray emission 
and tends to predic t spectral slopes which rise too s teeply 
in the UV ( see e.g. IAntonuccillT99l iBonning et "alll2007l : 
iDavis et alJl2007L and the AD model to SED comparison 
for the low M objects in Figure [7]). However, the model 
has had some success in other sources. It does a rather 
good job of reproducing the thermally dominant states 
of BH X-ray binaries, in which the emission is believed to 
come predominantly from an AD. Detailed spectral mod- 
els similar to the ones used here (jDavis fc Hubenvl[20?)fih 
not only reproduce the spectrum, but also the spectral 
evolu tion as M varies (e.g lDavis et al. 2006; S hafee et al.l 
[200fih . 

A central premise to our 77 estimates is that some mech- 
anism acts to redistribute flux from UV to the X-ray fre- 
quencies. In principle, such a mechanism can account 
for the mismatch between observation and theory in the 
UV while simultaneously explaining the larger than ex- 
pected X-ray flux. This could be related to advection, 
a Comptonizing coronae, a warm skin, etc. Regardless 
of the mechanism invoked, our rj estimates rely on the 
assumption that it does not act in regions where most of 
the optical photons are emitted. If some mechanism is 
acting on the optically emitting regions which breaks the 
assumption that local radiation flux (and work done by 
the accretion stress) balances the local release of binding 
energy, it could clearly modify our estimates. 

This standard disk model gives, at best, mixed results 
for explaining the sizes of continuum emission regions 
inferred from microlensing constraints in gravitationally 
lensed QSOs. The major discrepancy is that the opti- 
cal to UV emission is constrained to come from radii 
which are factors of ~ 3 — 10 larger than expected from a 



standard disk Ce.g. iMortonson et al.l 


20051 iPoolev et al.l 


120061 l2007HMorEan et al.l 20ia 120081: 


Dai et al.ll2010D. It 



also tends to give 77 estimates which are lower (log 77 
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— 1.77 + log than those derived from the Soltan argu- 
ment ijMorgan et al.|[2()10D . More optimistically, the rel- 
ative scaling of radius with observing wavelength (i? oc 
A*/^) is consiste nt with microlensing constraints in sev- 
eral s ources (e.g. lEigenbrod et al.ll2008HPoindexter et alj 
|2008() . Fu rthermore, the M dep endence of the microlens- 
ing radii ([Morgan et al.l 120101 ) are also consistent with 
theoretical expectations (R cx M^/^). 

In principle, one could reconcile the discrepancy in ab- 
solute microlensing radii by assuming a flux profile which 
falls off less steeply with radius. This would also improve 
the agreement between observations and spectral slopes. 
If we parametrize the radiative flux as oc R~^, the 
standard disk corresponds to /3 = 3 and the same argu- 
ments which leads to equation ([5|) yields 

OcM2/^M(2~4//^)i.(4-8//3) (18) 

As (3 decreases the spectral slope flattens, with i^L^, inde- 
pendent of v at frequencies that are unaffected by trun- 
cation at the inner or outer radii of the disk. Therefore, 
a model with /3 < 3 would give larger radii for the optical 
emitting regions and also flatter optical to UV spectral 
slopes, in better agreement with observations. 

Suggestively, a model with /3 closer to 2 would also give 
a ratio of iopt/^boi which is less dependent on AI at fixed 
rj. This can be inferred from equation ^TE\\ . which shows 
that the dependence of lyL^ on M weakens as /3 — >■ 2. For 
13 — 2, the M needed to match Lopt would be indepen- 
dent of AI and a constant rj model would be consistent 
with a constant Lopt/iboi- Such a scaling could, for ex- 
ample, result from irradiation. If the optical emission 
came from a flared region of the outer disk which was 
irradiated by the emission from the inner disk, the re- 
processed radiation could have a radial profile with /3 
closer to 2 at larger radii. Thus one might obtain a more 
constant rj by accounting for irradiation in the model if 
it dominates the intrinsic, local emission. 

However, there are two main problems with such an in- 
terpretation. First, it would give a scaling R oc A*/'^ that 
would not agree with microlensing cons traints, which as 
noted above are consistent w ith f3 ^ 3 (jEigenbrod et al.l 
120081: IPoindexter et al.l2008D . Secondly, if the sources are 
irradiated, we must be observing more optical emission 
than we would have from a bare AD. By using Lopt with- 
out accounting for the "extra" reprocessed emission we 
would have overestimated AI and underestimated rj. Al- 
though irradiation could then account for the low AI sys- 
tems that have low i], it could only increase our already 
uncomfortably high ry estimates for large M. Therefore, 
we do not believe that irradiation alone can explain the 
rj — AI correlation. 

If the assumed redistribution mechanism that shifts 
emission from the UV to X-ray depends mainly on r = 
R/Rg, it may be the case that the optical emission will be 
increasingly influenced by the mechanism as AI increases. 
This is plausible if the redistribution mechanism primar- 
ily acts in a confined region near the innermost radius of 
the disk, such as might be the case for Comptonization 
of disk photons in a compact corona. As AI increases 
the range of optically emitting radii moves to smaller r. 
Therefore a larger fraction of the optical emission may 
be subject to the redistribution, and the observed Lopt 



may be underestimating the true AI and overestimating 
rj. Such a scenario could potentially explain the high val- 
ues of rj at large AI, but cannot account for the low rj at 
small AI. 

Another possibility is that the assumption of a con- 
stant AI through the disk may be incorrect if there is 
outflow interior to the optical emitting region that car- 
ries away a sizable fraction of the accreting mass. Such 
an outflow converts some of the accretion power into 
mechanical luminosity, rather than radiative luminosity. 
The outfiow can be in the form of a wind, as seen in 
broad absorption line quasars, or as a jet, seen in radio 
loud quasars. A high mechanical luminosity can be ac- 
commodated in low M AGN, where the radiative rj is 
low. However, low M and low Lopt AGN generally do 
not sho w high velocity outfiow s, in particular in the PG 
sample (|Laor fc Brandt! [2002D . In high AI systems only 
a small fraction of the accretion power can go out as 
mechanical luminosity, as the radiative efficiency is close 
to maximal. However, both jets and fast UV absorbing 
outflows are commonly seen in high AI and high Lopt 
systems, suggesting they carry only a small fraction of 
the radiative power. Thus, mechanical outflows does not 
appear to go in the direction that can weaken the cor- 
relation of 77 with M. But, the high 77 values in high 
AI systems provides an interesting limit on the possible 
mechanical feedback of AGN on their environment. We 
note in passing that some of the AGN radiative and me- 
chanical luminosity may be obtained by tapping the BH 
spin, which can drive 77 to values above the theoretical 
AD upper limit. 

4.4. Evolution of rj with AI 

Thus far we have focused on potential errors in our as- 
sumptions which could artificially introduce correlations 
between rj and AI , but it is entirely possible that such a 
correlation is real. As noted previously, our mean log 77 
is consistent with estimates derived with the Soltan ar- 
gument and the range of inferred 77 is physically realiz- 
able. Given the uncertainties in AI and Lboi, it is entirely 
plausible that all 77 estimates are below the theoretical 
maximum for a rapidly accreting BH (a* = 0.998) with 
no torque on the inner boundary. If magnetohydrody- 
namic torques are pres ent, the spin may be limited to 
lower values (a* ~ 0.9, iGammie et al.l[200^ and refer- 
ences therein). However, torques at the inner edge of 
the disk cou ld, in principle, yield even higher efficiencies 
(|Agol k Kro lik 2000). 

If a* is the predominant factor which determines 77, 
as in the standard thin disk mo del (jShakura fc SunvaevI 
[T973I: INovikov fc Thornd[T973l) . one would expect 77 to 
vary with AI if a* varies with AI. In fact, detailed 
semi-analytic models of galaxy formation that simulate 
the accretion histories of super-mas sive BHs can, find 
a strong dependence of a* on M (iLagos et al.l [200I 
iFanidakis et al.ll2009D . For their Model A. ILagos et al.l 
(2009!) find a general trend of the average a* increasing 
with AI. For M < IO^Mq, they find a range of a^,, cover- 
ing the range from to 1, but with more systems at low 
a*. At M > IO^AIq, they find mo st systems with a^, ~ 1 . 
For their chaotic accretion model. iFanidakis et al.l ()2009l ) 
mostly find a* < 0.5 for M < 10^ Mq, with some systems 
counter rotating (o* < 0). At higher mass a* increases 
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Fig. 16.— Spins as a fun ction of M for the 80 PG QSOs in 
our sample. We use the the INovikov fc Thornel 1 11973) model to 
derive a* for each value of rj. Values of 77 > 0.31 or r; < 0.038 
are set to these limits, corresponding to a* = 0.998 and a» = 
— 1, respectively, and represented by filled triangles for plotting 
purposes. The massive BHs with logM > 8 generally have a, > 
0.5, and the lower mass BH s have < 0.5. For comparison we plot 
contours representing the IFanidakis et al.l |2009) a, distribution 
(see their figure 9). The thick dashed curve represents the median, 
while the shaded contours contain 60% and 80% of systems at the 
corresponding mass. This is broadly consistent with various earlier 
suggestions for a rise in a* with M . 



with a median a* ~ 0.6 - 0.8 fo r M > IQ- ^ MfD. The pro- 
longed accretion model of .Fanidakis et al.l (|2009[ ) yields a 
much different distribution, with a* ~ 1 for M < IO^Mq 
and dropping to a* ~ 0.8 — 0.9 at higher mass. 

In order to compare with these and other predictions, 
we convert 77 to a, for each source in our sample. This 
is done assuming disk model with no inner torque and 
an inner radius corresponding to rmg. The minimum 
and maximum rj correspond to 0.42 and 0.038 for mod- 
els with a* = 1 and -1, respectively. The results are 
plotted in figure 1161 In cases where our rj estimates lie 
outside the allowed range, we assume a, = —1 or 1 for 
r/ below and above the allowed range, respectively. The 
relatively gradual evolution in ry with M translates into a 
more abrupt evolution in a* due to the rapid rise in 77 as 
a* approaches unity. For comparison purposes, we have 
plotted dashed curves outlining the a* distribution of the 
IFanidakis et al.l (|2009t ) chaotic accretion model in figure 
1161 The thick curve is the median of the distribution and 
the shaded contours contain 60% and 80% of systems at 
the corresponding r nass. 

Both model A of iLagos et al.l (120091) an d the chaotic 
accretion model of IFanidakis et al.l ( 20091) are qualita- 
tively consistent with our results, in that the average a* 
increases with M. However, neither is a perfect match to 
our correlation. Our es timates wou l d not accommodate 
the high a* at low M in lLagos et al.l (|20Q9[ ). and our high 
M systems would r equire very near uni ty, so the spins 
in this range in the Fa nidakis et al.l (|2009[ ) model are not 
quite high enough. Of course, both of these calculations 
are subject to uncertainties in what they assume about 
the angular momentum of the accreting gas, so there may 
be some leeway to accom modate our resu l ts. T he pro- 
longed accretion model of IFanidakis et al.1 (|2009f ) would 
be completely inconsistent with our estimates. 

Further support for spin evolution comes from mod- 



els for radio emission in Q SOs that require high a* (e.g 
IBlandford fc ZnaieklUlTTl) . Since radio loud QSOs are 
associated with large M, high a* would be required at 
large M if these models are correct. The minimum value 
of a* needed for significant jet production is somewhat 
uncertain in these models, so it is difficult to make this a 
quantitative argument, but it is qualitatively consistent 
with our finding of high 77 at large M. 

Other factors besides a* may also play a role in set- 
ting rj. One possibility is that advection may be op- 
erating near the inner-m ost radii for system s accreting 
near the Eddington rate (|Abramowicz eTaD 11988). Fig- 
ure [T3| shows our Lboi/^Edd distribution as a function 
of M. Most of the systems with iboi/^Edd > 1 have 
M < lO^A^Q. This may, in part, be due to the low- 
est M systems having their M predominantly underes- 
timated as discussed above. Alternatively, it may be a 
selection effect resulting from the requirement that the 
QSO outshine its host galaxy, which will require a higher 
^boi/^Edd in disk dominated galaxies. If these systems 
are truly accreting at or above the Eddington rate, the 
advection time is comparable to the photon diffusion 
time and a substantial fraction of the radiation might get 
advected into the BH before radiating. If these sources 
are accreting at only mildly super-Eddington rates, ad- 
vection should have no effect at larger radii where the 
optical is emitted, and our M estimates should still be 
reliable. These systems would then be truly radiatively 
inefficient. Comparison of figures |9| and [TSj show that 
the effect of advection in flows with Lboi/-^Edd ^ 1 could 
plausibly account for some, but not all, of the lowest 77 
sources. 

5. CONCLUSIONS 

We have estimated the accretion rate by fitting radia- 
tively efficient AD model SEDs (hereafter the standard 
model) to the optical emission in a sample of 80 PG 
QSOs. This method is insensitive to properties of the 
accretion fiow and BH spacetime near the inner edge of 
the disk (torques at the inner radius, BH spin, advection, 
etc.) as long as the emission comes from large radius. We 
use detailed AD model SEDs which are computed from 
non-LTE atmospheres and include relativistic effects on 
photon geodesies, but find our results are qualitatively 
reproduced by simple, non-relativistic local blackbody 
relations. The derived accretion rates are nearly insen- 
sitive to spin at low black hole mass (Af < IO^A/q) and 
only weakly sensitive at higher mass when the fitting 
is done at optical frequencies. The accretion rates are 
more sensitive to assumed BH masses, which are esti- 
mated using broad line region virial methods as well as 
masses derived from the A/ — relation for 13 sources 
with bulge velocity dispersion measurements. 

Our sample of 80 PG QSOs was chosen because they 
have ample broadband coverage at optical to far UV 
and X-ray frequencies. This allowed us to robustly esti- 
mate the bolometric luminosity modulo some uncertainty 
in the unobservable EUV emission. These luminosities, 
combined with our estimates of the accretion rate, allow 
us to compute the radiative efficiency for each source in 
our sample. We find a mean efficiency of log 77 = — 1±0.5, 
in agreement with integral constraints derived by match- 
ing local black hole mass density to the integrated quasar 
luminosity function (i.e. the Soltan argument). This ba- 
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sic agreement suggests that the standard model can pro- 
vide a reasonable first-order approximation to real ac- 
cretion flows, at least at radii where most of the optical 
emission is produced. 

We find a strong correlation of efficiency with BH mass 
(approximately 77 oc M^/^ in our sample) extending from 
77 ~ 0.01 for low masses to to 77 < 1 at the highest 
masses. This relation arises because the ratio of the op- 
tical to bolometric luminosity is roughly independent of 
BH mass, whereas a constant efficiency thin disk model 
would predict a substantial increase in this ratio as mass 
increases. 

We consider three possibilities for explaining the 
77 — M correlation; 

1) The correlation is real. It could plausibly arise 
from a mass dependence of the BH spin driven by the 
differing accretion histories of black holes at different 
masses. Semi-analytic models of galaxy formations that 
attempt t o model spin distributions of su permassive 
BHs (e.g. iLagos et al.l [20091: iFanidakis et all boOQ) find 
spin dependencies which are qualitatively, although not 
quantitatively, consistent with our trend of increasing 
efficiency with increasing BH mass. 

2) One or more of the observables or input parameters 
to the model are incorrectly estimated. Indeed, scatter 
in the broad line region based BH mass estimates 
probably contributes to some of the trend. We argue 
that the highest and lowest BH masses will tend to be 
overestimated and underestimated, respectively. This, 
in turn, leads to overestimates of the efficiency at the 
highest masses and underestimates at lowest masses. 
However, it would be difficult for these correlated 
errors to explain all of the trend, unless the true mass 
distribution is very narrow (< 1 dex), which seem highly 
unlikely given the large range of bolometric luminosities 
in our sample and the significantly larger range of M 
derived in other studies of AGN. 

3) The standard (bare) AD model does not adequately 
approximate the dependence of real accretion flows 
on mass. It has well-known difficulties reproducing 
the observed SEDs of AGN and microlensing sizes of 
emission regions. A central assumption of this work is 
that standard model works well at larger radii where 
the optical emission is predominantly produced, but 
that some mechanism operates very near the black hole 
to redistribute flux from UV to X-ray frequencies. It 



is possible that the standard model fails in the optical 
emitting regions as well and that real flows naturally 
give rise to a relatively constant ratio of optical to 
bolometric luminosity which does not depend strongly 
on mass. We consider a number of modifications to 
the standard model, including irradiation. However, 
no single possibility we considered seemed capable of 
explaining the whole trend at both low and high masses, 
while maintaining a nearly constant 77 whose value was 
consistent with Soltan argument. 

The role of errors in the mass estimates could be defini- 
tively addressed by more precise mass estimates. This 
would either require reducing the scatter in the BLR 
estimates or obtaining a larger sample of sources with 
well-measured bulge velocity dispersions. Differentiating 
the degree to which the trend reflects real evolution in 
efficiency with mass or inapplicability of the standard 
models assumptions is more difficult. Microlensing mod- 
els can provid e independent estim ates for the efficiency in 
lensed QSOs ([Morgan et al.ll20fo( ). although only a hand- 
ful are currently available. Alternatively, relativistically 
broadened Fe Ka lines could provide spin estimates from 
which we could infer efficiencies, but (again) precise con- 
straints are only available for a few sources. Some light 
on what controls rj may be shed by the variability of 7/ 
on timescales longer than the viscous timescales, when 
the AD may be quasistatic. If Lboi oc L^p^, then 77 is 
constant, which is consistent with being spin driven. Un- 
til such studies are available for a significant sample of 
sources, this will likely remain an open question. 
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TABLE 1 — Continued 



Object 


H/3 FWHM^ 


logLopt'^ 




log M'^ 




log Lboi'' 


logr? 


2214 + 139 
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-0.50 
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-1.01 


2251 + 113 
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9.43 
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-0.27 



Note. — Summary of M and r] derived using Mblr estimates. We report M and M without error as systematic uncertainties in 
the estimation methods dominate the statistical uncertainty in the input data. These uncertainties are discussed further in the text. For 
brevity, we do not report the uncertainties in logr; because they are identical to the uncertainties in logLboi. 

The FWHM of the H/3 line in units of km s"!. 
^ Lopt and Lbol measured in units of erg cm~^ s~^. 
^ M measured in units of Mq. 

M measured in units of Mq yr~^. 

TABLE 2 



Object log Mo-" logAf*^ logj? 



0007 


+106 


8.14 


-0.26 


-0.88 + 0.17 


0050 


+124 


8.02 


-0.36 


-1.18 ±0.04 


1119 


+120 


7.76 


-0.73 


-0.74 ±0.34 


1126 


-041 


8.08 


-0.75 


-0.75 ±0.28 


1229 


+204 


7.76 


-0.38 


-1.22 ±0.04 


1302 


-102 


9.09 


0.62 


-0.76 ± 0.12 


1309 


+355 


8.42 


0.08 


-1.10 ±0.05 


1404 


+226 


8.43 


-1.07 


-0.38 ±0.26 


1426 


+015 


7.99 


0.11 


-0.92 ± 0.24 


1617 


+175 


7.97 


0.01 


-1.23 ±0.14 


2130 


+099 


7.87 


-0.31 


-0.82 ± 0.32 


2214 


+139 


7.70 


-0.14 


-1.37 ±0.08 


1444 


+407 


8.71 


0.06 


-0.44 ±0.28 



Note. — Summary of M and derived using Mo- estimates. We report M and M without error as systematic uncertainties in the 
estimation methods dominate the statistical uncertainty in the input data. These uncertainties are discussed further in the text. 
^ M measured in units of Mq. 
^ M measured in units of Mq yr~^ . 



